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Abstract 

This paper studies the hnk between the number of critical eigenvalues and the number of 
delays in certain classes of delay-differential equations. There are two main results. The first 
states that for k purely imaginary numbers which are linearly independent over the rationals, 
there exists a scalar delay-differential equation depending on k fixed delays whose spectrum 
contains those k purely imaginary numbers. The second result is a generalization of the first 
result for delay-differential equations which admit a characteristic equation consisting of a 
product of s factors of scalar type. In the second result, the k eigenvalues can be distributed 
amongst the different factors. Since the characteristic equation of scalar equations contain only 
exponential terms, the proof exploits a toroidal structure which comes from the arguments of 
the exponential terms in the characteristic equation. Our second result is applied to delay 
coupled D„-symmetric cell systems with one-dimensional cells. In particular, we provide a 
general characterization of delay coupled D„-symmetric systems with arbitrary number of 
delays and cell dimension. 



1 Introduction and Background 

Delay-differential equations (DDEs) have been used as mathematical models for phenomena in 
population dynamics [20] , physiology [12l [3] , physics [22] , climate modelling |28j and engineer- 
ing [28] amongst others. Delay-differential equations behave like abstract ordinary differential 
equations (DDEs) on an infinite-dimensional (Banach) phase space and many results which are 
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known for ODEs on finite-dimensional spaces have analogs in the context of DDEs. Many scalar 
delay-differential equation models have been developed over the years such as for Cheynes-Stokes 
respiration [12] and the regulation of hematopoiesis [12], the delayed Nicholson blowflies equa- 
tion [17] in population dynamics, a two-delay model of an experiment on Parkinsonian tremor [2] 
and many more. 

The bifurcation analysis of DDEs is done essentially in the same way as in ODEs, although 
the technical details differ. Consider the neighborhood of an equilibrium solution of a nonlinear 
DDE, then the analysis of the linearization at the equilibrium point leads to stable, unstable and 
centre invariant subspaces where only the stable subspace is infinite-dimensional. There exists 
local invariant manifolds (stable, unstable and center manifolds) tangent to the corresponding 
invariant subspaces of the linearized equations about the equilibrium point on which the flow 
near the equilibrium is either exponentially attracting (stable manifold), exponentially repelling 
(unstable manifold), or non- hyperbolic (center manifold). Now, bifurcations near equilibria are 
determined by the flow on the centre manifold and the dimension of this manifold is determined 
by the number of eigenvalues of the linearization on the imaginary axis. 

The first result of our paper is Theorem 12.11 and goes as follows. Consider n nonzero imagi- 
nary numbers iivi, . . . ,iuJn where the imaginary parts uji, . . . ,ujn are positive and not rationally 
dependent. We show that there exists a scalar linear delay-differential equation depending on n 
discrete delays written 

n 

X = ^ajx{t-Tj) (1) 

i=i 

where x G R, Oj E M and Tj G [0,t] for all j = l,...,n such that the characteristic equation 
of (dl), given by 

n 

A-^a,e^-^=0, (2) 
i=i 

has eigenvalues ztiwi, . . . , ±iujn- This result generalizes explicit computations done in the case of 
one and two delays, see [191 El I]- The proof is done by embedding the problem as a mapping 
which is solved by the implicit function theorem at a carefully chosen point. From the implicit 
function theorem, we are able to define a smooth mapping whose transversal intersection with 
a dense curve on an n-dimensional torus provide solutions. The incommensurability of the n 
frequencies enables us to define the dense curve on the n-torus. This type of argument using a 
dense curve on an n-dimensional torus was used in Choi and LeBlanc [7|. 

This result falls within the category of so-called realization theorems. For instance, the re- 
alization theorem of linear ODEs by linear DDEs obtained by Faria and Magalhaes [llj. They 
show that for any finite dimensional matrix B, a necessary and sufficient condition for the exis- 
tence of a bounded linear operator Co from C([— r, 0], M") into M" with infinitesimal generator 
having spectrum containing the spectrum of B is that n be larger than or equal to the largest 
number of Jordan blocks associated with each eigenvalue of B. Other results in this direction are 
concerned with the realization of finite jets of ODEs on a finite-dimensional centre manifold by 
delay-differential equations, see |1H 17]. To our knowledge, the realization theorems in this paper 
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are the first general results linking the number of critical eigenvalues of linear delay-differential 
equations with the number of discrete delays. 

The next significant result is an openness theorem, that is, the realization of n imaginary 
numbers (not necessarily rationally independent) as eigenvalues of a linear scalar delay-differential 
equation is valid in a neighborhood of any set of n rationally independent imaginary numbers. 
The proof of this theorem also relies on the implicit function theorem. 

We then turn our attention to the context of symmetric systems of delay-differential equations. 
Several examples of symmetric systems of DDEs \15\ [23] have characteristic equations which 
decompose in factors, some of which have the same form as the characteristic equation ([2]). The 
decomposition of the characteristic equation is induced by the isotypic decomposition of the space 
and we present a general derivation of this decomposition. We show that isotypic components 
consisting of a unique one-dimensional complex irreducible representation contribute a factor of 
the form ([2]) in the characteristic equation and so Theorem 12.11 can be applied directly to each 
such factors separately. 

We present a generalization of Theorem 12.11 to the case where several factors of the charac- 
teristic equation have purely imaginary eigenvalues simultaneously. Theorem 12.41 shows that a 
set of n rationally independent purely imaginary complex numbers can be realized from several 
factors of the characteristic equation of a delay-differential equation with n delays given some 
nondegeneracy conditions on the characteristic equation are satisfied. The statement of the the- 
orem is independent of any symmetric structure and the proof is a generalization of the proof of 
Theorem 12.11 

We illustrate the above result on D„-symmetric rings of n delay equations with delayed 
coupling. Hopf bifurcation from such symmetric networks have been studied by several au- 
thors [6l [Ml [m [TU [211 ESI |29j. In order to apply Theorem 12.41 to this context, we derive an 
explicit form of the coupling matrix in terms of the connections in the graph representation of 
the ring for cells of any dimension and arbitrary numbers of connections and delays. This is a 
generalization of the networks considered in the articles listed above in this paragraph. 

We specialize to the case of one-dimensional cells and we can then obtain general formulae 
for the factors of the form ([2]) in the characteristic equation of the D„-symmetric ring. Given 
a nondegeneracy assumption of Theorem 12.41 is satisfied, then it can be applied directly to D„ 
symmetric coupled cell systems with n odd. We do not treat the case n even because the form of 
the equations needed for Theorem 12.41 is not satisfied. We illustrate this fact in a D4-symmetric 
example. 

The paper is organized as follows. The first section contains brief preliminary remarks and then 
we state and prove our main result (Theorem 12. ip and the openness result (Theorem 12. 2p . Then 
we introduce the context leading to Theorem 12.41 and state this result. Section [3] is devoted to F- 
symmetric systems of delay-differential equations and the section begins with a general discussion. 
Section 13.11 presents a characterization of D„-symmetric rings of delay coupled cells with an 
arbitrary number of delays and derive the characteristic equation in the case of one-dimensional 
cells. In Section [3.21 Theorem 12.41 is applied to D„-symmetric rings of one-dimensional cells with 
n odd. Section [H has the proof of Theorem 12. 4[ We conclude by a discussion of open problems 
along the lines of the ones presented in this paper. 
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2 Realization theorems 

We now discuss some aspects of the spectral theory of hnear scalar delay-differential equations. 
In fact, we just introduce the basic facts, in a non-abstract setting, needed for the statement of 
our first main theorem. For a complete treatment, see Diekmann et al [9j or Hale et Lunel [19] . 
Consider the scalar delay-differential equation 

n 

xit) = ^ajx{t-Tj) (3) 

i=i 

where aj G M and tj £ [0,r] for all j = 1, . . . ,n and r > 0. The characteristic equation for ([3]) 
can be obtained by substituting x{t) = Ce'^^, where C is a constant, into the equation. Thus, 

and by rearranging the terms we obtain 

\ 

3=1 J 

So, x(t) is a nonzero solution of ([3]) if and only if 

n 

A(A) ■= X-^^aje-^^^ =0. 
i=i 

The complex number A is an eigenvalue of equation ^ if it is a solution of the characteristic 
equation A (A) = 0. 

The question we address in this paper is related to the number of imaginary eigenvalues (with 
incommensurable frequencies) which can satisfy A (A) = 0. The case n = 1 with one nonzero 
delay is a straightforward calculation and A (A) = for only one nonzero imaginary eigenvalue 
A, see [H]. The case n = 2 with ri = and T2 S (0,t] in ([3]) can be found in [9]. There, it is 
shown that A(A) = can have at most two nonzero imaginary eigenvalues. The case n = 2 with 
Ti,T2 > is done in [1] where it is shown that A(A) = can have at most two nonzero imaginary 
eigenvalues. We are now ready to state our first result. 

Theorem 2.1 Suppose wi > 0, 0^2 > 0, . . . , cj„ > are linearly independent over the rationals. 
Then there exists ri > 0, r2 > 0, . . . , r,i > 0, ai € M, 02 E M, . . . , G K such that the linear 
delay differential equation 

x{t) = ai x{t - Ti) + a2x{t - T2) -\ h a„ x{t - r„) (4) 

has solutions x^it) = e^*^-'* for all j = 1, . . . , n. 
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Proof: A necessary and sufficient condition for tlie conclusion of tfie tlieorem to liold is tliat 
the following algebraic system of 2n equations 



E 

k=l 
n 

k=l 



UjJ 



J 



1, 



-lU) 



31 



1, 



(5) 



has a solution in the 2n unknowns (n, r2, . . . , Tn, ai, 02, . . . , an). Although ([5]) is in complex form, 
since the second equation in ([5]) is just the complex conjugate of the first equation in ([5]), system 
([5]) is equivalent to a system of 2n real equations. This fact is taken for granted throughout the 
sequel, even though we continue to use complex notation. 
It is useful to use the following matrix notation for ([5]) 



P(t;u;) 
P(-t:uj) 



(6) 



where, uj = (wi, . . . , A = (ai,...,an), superscript T denotes transpose, and P(r;a;) 
P(ri, . . . , r„; wi, . . . , LOn) is the n x n matrix whose entry at row j column k is 



:P(r;a;)] 



jk 



Note that P{t]lo) = P{-t;uj). 



Remark: Let V denote the 2n x n matrix of coefficients of the left-hand side of // we view 
Ti G M'^ , . . . ,Tn € M"*" as free parameters for these coefficients, we see that each column of the 
matrix V is dense on an n-torus, T = (S"*^)", since the iOj are rationally independent. This remark 
is a key point in the following discussion. 



Define V = T"' = ((S^)"')'^, and choose coordinates on V as follows: 

V = {cI> = (cDi,...,$«)|c^^ = (<pj,... VjG(S^r,J = l,...,n}. 
Consider the following mapping associated to ([6]): 

F : V X M" ^ 

defined by 



t>2n 



P{<f) 
Pi-<l>) 



A^-i 



ijj 



-UJ 



(7) 



where A and uj are as previously defined, and -P(<^) is the n x n matrix whose entry at row j 
column k is 



P{^) 



p(q>^ 



jk 



e . 
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For reasons that are clear from the sequel, we choose to single out the last torus in the product 
V = T", and use the variable ^ as coordinates for this torus, i.e. 

^ = (^l,...,Vn) = 1>" = (<^?,...,V'^). 

Thus, we write V = V$ x Vij-, where = T""-*^ and Viji = T, so that 

F : V$ X VvK X M" I — > M^n 

is written as F(^, A;uj) in d?]) (we have re-labeled <I> = (<I>^, . . . , (J)""^) to designate coordinates 
for V$ ^ T"-i). 

Definition: If {ei, . . . ,6^} denotes the canonical basis of vectors in M", we define the vectors 

n 

vi, . . . ,Vnhy vi = Ck, and for j = 2, . . . , n, 

k=l 

e=o 

By construction, the set {vi, . . . ,Vn} is linearly independent, and so the n x n matrix I, whose 
j^^^ column is the vector vj, is invertible. We also define the n x n invertible matrix Uj to be the 

diagonal matrix whose diagonal element is the k^^ component of the vector vj (in particular, 
Ui is the identity matrix). Note also that U^^ = Uj, j = 1, . . . ,n. 



Consider the following point in V$ x Y^,: 
then it is easy to compute that 

P($,$)=iJ, (8) 

where P is as in ([7]) . If we define 

= {ai,...,anf = I-^LO^ 

then 

F{$,^,A;lo) = 0. 

Because the ivj are rationally independent, it follows that the components ak of A are all non-zero. 
Furthermore, it is easy to compute the following derivatives 

^ ( CLnUn il \ 

J = D^^^A-^F{^,^,A-oj)=[ 

\ anUn -iT / 
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and 

^ ^ ^ ( aiUl a2U2 ■ ■ ■ Cln-lUn-l 

K = D^F{^,^!,A-Lo) = 
The matrix J is invertible, since its inverse is easily computed as 



J- 



1 7/ 1 7/ 

2 2 



By the impHcit function theorem, there exists a neighborhood of <I> in V$ and a unique smooth 
function 

G : iV I — > Vvt X M" 

1 — > G($) = (G^($),Ga($)) 

such that 

G($) = ($,l) 

and 

F($,G($);l<j) =0, V«> G TV. (9) 
Imphcit differentiation of ([9]) yields that 



DG(<^) = ( I = -J-^K 

DGa{^) 



(10) 



V ••• 

where denotes the n x n zero matrix. Consequently, 

DG^i^) = -I^UnUi -^UnU2 ■■■ -^UnUn-l ), (H) 

and it follows that the mapping 

G^:N — ^ V* 

is regular at ^. 

Consider the following n — 1 vectors in (M" )"~"^: 

Wi = (a;,0,...,0,0) 
W2 = (0,w,0,...,0,0) 



Wn-2 = 
Wn-l = 



(0,0,...,w,0) 
(0,0,..., 0,u;), 
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where represents the vector in M", and we recall that lo = {ui, ... , a;„). 
The set {Wi, ... , W„_i} is linearly independent, so for any $ G V$, the set 

ji-i 

5$ = { $ + ^ (Tj Wj (mod V$) I < \aj\ « 1, j = l,...,n} 

is a small n — 1-dimensional surface through $ in V$. We are interested in showing that for $ 
close enough to $ in A'^, the image of <S$ by in V* is transverse to the integral curves of the 
vector field * = a;. To show this, we consider the function 

T:N — >R 

defined by 

T($) = det {DG^{<^) • W[ DG^{<^) -W^ ■■■ DG^{^) ■ W^_^ a;^ ) . (12) 
Obviously, T is continuous, and 

r($) = det( -ll-UnUiiO^ -^UnU2U;^ ■■■ -^UnUn-luF UnUnUJ^ ) 

= T — detW„det( aiUiuF 0,21(200^ ••• hn-iUn-ioj^ UnOJ^ ) 

h 

(wi UJ2-- - iOn) (ai 0.2 ■■■ a„_i) 
= detX 

7^ 0. 

It follows that there is a neighborhood N' C N in which T 7^ 0. This is equivalent to saying that 
for all $ G iV', the image of 5$ by G^ in is transverse to the integral curves of the vector 
field ^' = w. 

For each j = l,...,n — 1, the integral curves of the vector field = ui are dense. Thus, for 
any e > 0, there is a Tj^^ > and an Sj^s > such that the integral curve segment 

{ = TjUJ (mod (S^)") I Tj, , - Sj,, < Tj < Tj,e + Sj,e } 

is in the e-ball centered on ^vj in (S^)". For e > small enough, the surface 

{ $ = (riW, T2UJ, Tn-ILO) (mod V$) I Tj^e - Sj,e < Tj < Tj^e + Sj,£ } 

is contained in A^' and coincides with the surface for 

^* = (ti.eO;, T2,£a;, . . . , T^n-i),e^) (mod V$). 

Thus, by our previous result, the n — 1-dimensional surface G^(»S$*) is transverse to the integral 
curves of ^ = w in Y^,. Since these integral curves are dense in V^, there are infinitely many 
intersections with G^((S$*) near the point ^ = G^($). 
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Let G^{S^*) be such an intersection point near ^. Then there is a r„> such that 

° o 

vl/=r„ io (mod Vij,). 

o o o oo o 

Let $€ 5$* be such that Gii,{^) =^1/. Then there are ti> 0, T2> 0, . . . ,t„_i> such that 

o o o o 

(j)= (ti r2 . . . , r„_i w) (mod V$). 



It follows from dH) that F(<I), xj/, G'a(<^); = 0, and by construction, this corresponds to a 
solution of dll). ■ 

The next theorem shows that the previous realization result holds for open sets near solutions 
found in Theorem 12.11 

Theorem 2.2 Suppose uji > 0, L02 > 0, . . . ,ujn > are linearly independent over the rationals. 
There exists a neigborhood M of u) = (wi, . . . , in M" and a smooth mapping 



H : V 



UJ 



H{uj) = {t{uj),A{uj)) = {{ti{uj),. . . ,Tn{uj)),{ai{u;),. . . ,an{uj))) 



such that 



k=l 
n 

k=l 



,n 



iUJj, j = 1, 

-iujj, j = l,...,n 



(13) 



for all UJ £ J\f . 



Proof: We consider the system F = given by ([7]) . We have already shown in Theorem 12.11 
that, for fixed uj linearly independent over the rationals, there exists infinitely many solutions 
to -F = 0. We again use an implicit function theorem argument combined with the density of 
irrational torus flows. 
Consider the mapping 

Q : M'" X M'" X R" — > M" x 

(14) 

{t,A,L0) I > Q{T,A,Lo)=F{{TlUJ,...,Tn-lL0),TnL0,A]L0), 

where F is as in ([7]). Therefore 

/ • 



DtQ{t, a, uj) = L'(($i^...^$n-i)^^)F((ria;, . . . , Tn-iuj) , TnUJ , A; 



UJ 



uj^ 



\ 




\ ■■■ uj^ J 



(15) 
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where each in the matrix above is an n-dimensional zero column vector; and 

DaQ{t,A,u;) = DaF{{tiu;, . . . , T„_ia;), t„w, A; w). (16) 

Thus, we wish to show that the 2n x 2n matrix 

{DrQ{T,A,Lo) DaQ{t,A,lo)) (17) 

is invertible at the solutions to ([5]) we have found in Theorem 12.11 

For positive integers p and q, let Matp^g denote the space of p x q matrices. Consider the 
following mappings associated to (fTSi) . (fT6]l and (fT7|) : 



7^l : V$ X Vvt X M" X 



Mat2n,n 



defined by 



/ cj^ 



tj'^ 



7^l($, ^, A, Lo) = L)($,v,,)F($, ^, A; 



defined by 

and 

defined by 



\ 

7^2 : Vvt X V* X R" X R" — > Matsn.n 

n2{^,^,A,uj) = DaF{^,^,A;uj), 
7^ : V$ X X R" X R" — > Mat2n 2n 



\ 




n{^,^,A,uj) = {ni{<p,^,A,uj) n2{'p,^,A,uj)). 

Now, a simple computation (similar to those done in the proof of Theorem 12. ip shows that 



Z il 
Z -iX 



where 



If none of the aj vanish, then the n x n matrix Z is invertible, since its determinant is 



det Z = ajUJj detJ 7^ 0. 
i=i 
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Thus, 

vr 

By continuity, there is thus a neighborhood Af of the point — — {vi,..., v„) in V$ x V^r in which TZ 

is invertible. By Theorem 12 .11 there are infinitely many solutions of Q = (see (|14p ) in M, and the 
Jacobian matrix (I17p is thus invertible at these solutions. We get the conclusion of Theorem 12.21 
by the implicit function theorem. ■ 



2.1 Example: Da-symmetric system 

Theorem 12.11 is written in the context of scalar delay-differential equations. However, in this 
section, we look at an example of a Ds-symmetric system of delay-differential equations where 
Theorem 12 . 1 1 can be applied and then proceed to explain the generalization of this theorem which 
has applications to symmetric systems of delay-differential equations. 

Example 2.3 Let F = D3 be the group generated by k and 7 act on as follows: 

K.{xi,X2,Xs) = {xi,X3, X2), 7.(xi, X2, X3) = (X3, Xi, X2). 

Consider a linear Ds-symmetric coupled cell system with delayed coupling where each cell is 
one-dimensional and has the following form. 

Xi = aiXi{t - Ti) + a2[x2it - T2) + X3{t - T2)] 

X2 = aiX2it - Ti) + a2[x3{t - T2) + Xi{t - T2)] (18) 
X3 = aiX3{t - Ti) + a2[xi{t - T2) + X2{t - T2)]. 



where G M for i = 1,2,3 and 01,02,03 € M. The characteristic equation of system (jlSp is 
obtained by substituting {xi,X2,X3) = {wie^^ ,W2e'^^ ,W3e^^) into the equations. We obtain after 
simplification 

\wi = aie~'^'^^wi + a2e~^'^^ [w2 + W3] 
\w2 = aie~'^'^'^W2 + a2e~^'^'^ [W3 + wi] 
Xw3 = aie~'^'^^W3 + a2e~^'^^ [wi + W2] 

and rearranging the terms we have 

/O 1 l\"|/ii;i\ 

= 0. (19) 

where / is the 3x3 identity matrix. Letting a = \ — oie"'^'^^ and j3 = —026"'*''^^ equation p9p 
becomes 

/a f3 (3 \ / wi \ 

0. 





^0 


1 




/ Wi 


(A-oie-^^i)I-02e-^"'2 







:)] 


W2 






1 




\ 'W3 



a 


/3 




1 Wl 




a 




W2 






a 
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Let 

I a (5 f3\ 
A(A) =\ (3 a 13 \. 
\ p (3 a J 

We complexify and look at the isotypic decomposition of C'^ by the action of D3: 

C^ = Vo(BVi(B V2 

where Vq is the trivial representation of D3 and Vi , V2 are the standard irreducible representations 
of D3 (all representations are one-dimensional complex). A basis for Vq is uq = {v,v,vY, a basis 
for Vi is ui = {v,e^'^^^^v,e'^'^^^^v) and a basis for V2 is U2 = {v , e'^'^^^^v , e^^^^^v) . Therefore, 

A(A)'Uo = (a + 2p)uo 

and 

A(A)tii = (a — P)ui, A(A)ii2 = (a — (3)u2 

since e^'^*/^ = e^'^*/^. Therefore, in the basis given by the isotypic decomposition of C^, A (A) 
block diagonalizes so that we have 

a + 2(3 \ 

a-p ^2 = 0. 

a- p / \ ws / 

Hence, the eigenvalues are solutions to 

det A(A) = (a + 2p){a - pf = (A - aie'^^' - 2a2e-^^^){\ - aie'^^' + a2e-^^^f = 0. 

Each factor of the characteristic equation is of the same form as the characteristic equation for a 
scalar delay-differential equation. Therefore, by letting oi = oi and 02 = 2a2 in (A — aie^^'^^ — 
2026"^"^^), Theorem 12.11 applies directly. The same is true for the factor (A — aie~^'^^ + 026^'^'^^) 
where we let oi = oi and 0.2 = —02- Hence, for any choice of a set of complex numbers A = 
{iLVi,iL02} with 0^1, 11^2 > and rationally independent, there exists a linear D3 symmetric coupled 
cell system including A in its spectrum. 

In the context of bifurcation theory, the symmetry properties of the critical eigenspace depends 
on which factor contains the critical eigenvalue and this leads to different bifurcation behaviour. 
Two imaginary eigenvalues in the first factor corresponds to a nonresonant Hopf/Hopf mode 
interaction (without symmetry) while the second case leads to a nonresonant D3 Hopf/Hopf 
mode interaction. Details of the unfolding of these bifurcations can be found respectively in 
Kuznetsov [21] and Golubitsky et. al. [T3j. 

Note that Theorem 12.11 is not sufficient to guarantee the existence of a linear D3 symmetric 
coupled cell system with iuji satisfying the first factor and iijj2 satisfying the second factor simul- 
taneously. We characterize this situation as follows. Let h\ = h\ = 1 and h\ = 2 and 63 = ~1 ^-iid 
for fixed rationally independent iuJi,iuj2 (with uji,uj2 > 0), we look for ai,a2 and ti,T2 such that 
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and their complex conjugate equations are satisfied. This is the context of the next theorem which 
is a generaUzation of Theorem 12. 11 We state this result in a general form below and postpone the 
proof to Section 2] as it follows similar steps as the proof of Theorem 12.11 

Note that in the proof of Theorem 12.11 the matrix X defined in ([8]) is nonsingular by construc- 
tion and this is a crucial step in the argument. For this more general result we shall present, the 
matrix which holds a similar role is denoted by Zb since it is a matrix consisting of it the constants 
iPj^ which appear in equations (I20p . The form of this matrix is not relevant for the moment and 
the structure of the matrix is described in Section HI We are now ready to state the theorem. 

Theorem 2.4 Consider the factors 

iiU-ji^^K^'A (21) 
j=i \ k=i / 

of a characteristic polynomial where the constants 5-^ G R \ {0} are fixed for all j = 1, . . . ,r, 
k = 1, . . . ,n and suppose that detT^ / 0. Suppose that > 0, u;2 > 0, . . . , u;„ > are linearly 
independent over the rationals. Then there exists ti > 0, T2 > 0, . . . , r„ > 0, ai G M, 02 G 
M, . . . , On G M such that for all j = 1, . . . ,r, 

(^^-E^kKe'^""^ =0 

has roots iujj for (. = 1, . . . ,lj where li + ■ ■ ■ + ir = n. 

This theorem is applied in the following sections to the case of D„ symmetric coupled one- 
dimensional cell systems. If n odd, it is easy to show that 7^ holds, but for n even, several of 
the 6^^'s can be zero and Theorem 12.41 cannot be applied directly. 



3 Linear r-symmetric delay-differential equations 

For the results of this section, we find it convenient to introduce the well know abstract setting, 
see for instance Hale and Lunel [19], adapted to the symmetric case. Let C„ = C([— r, 0],C") be 
the Banach space of continuous functions from the interval [— r, 0], into C" (r > 0) endowed with 
the norm of uniform convergence. Consider the linear homogeneous RFDE 

m = Coizt), (22) 

where Cq is a bounded linear operator from C„ into C". We write 

AM = r drjiOMO), 
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llr ldr]{e)e 



e 



7 



where 77 is an n x n matrix- valued function of bounded variation defined on [— r, 0]. The charac- 
teristic equation is 

^0 

detA(A) = 0, where A(A) = A/„-y dr]{9)e^^ , (23) 

where I„ is the n x n identity matrix. Note that e^^ = e^^In- 

Suppose that F is a compact group of transformations acting Unearly on C". We say that (f22]l 
is T-equivariant if 

7 • r]{d) = r]{e) -7, V7 G r, e [-r, 0]. (24) 
The group action of F on C" induces an isotypic decomposition of C"": 

c" = e e • • • Vfc 

where Vi = Ui (B ■ ■ ■ Ui for irreducible representations Ui of T and Ui 9^ Uj for i ^ j. Since r}{6) 
commutes with the action of F, then 

v{0)Vi C Vi 

for all i = 1, . . . ,k. 

Therefore, A(A) also commutes with the representation of F. Indeed, for all 7 € F 

A(A)7 = A/7 
= 7A/ 
= 7A/ 

7(AI-7°^dr?(e)e^^)" = 7A(A). 
Thus, 

A{X)Vi C Vi 

and for all i = 1, . . . , /c and we can write A (A) in block diagonal form: 

A(A) =diag(Ai(A),...,Afc(A)). 
The characteristic equation then becomes 

k 

detA(A) = JJdetAi(A). 

i=l 

Therefore we are led to the following result. 

Proposition 3.1 Suppose that Vi = Ui and Ui is a one- dimensional irreducible representation of 
F . Then 

e 

detAi(A) = A-^Oj-e^^^- 
i=i 

Corollary 3.2 Theorem \2.1\ applies to factors of the characteristic equation which correspond to 
the context of Proposition \3.1\ 
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3.1 Delayed coupled cell systems with D„-symmetry 

Multiple authors have studied Hopf bifurcation in D„ symmetric rings of cells with delayed cou- 
pling where each cell is one-dimensional. The differential equation systems in those papers have 
the following general form. For i = 1, . . . ,n, the dynamics of cell i is given respectively for n odd 
and n even by: 

Xi{t) = f{Xi) + g{xi+i, Xj+(„_i)/2, Xi_(„_i)/2, . . . , Xi-i) (25) 

Xi{t) = f{Xi) +g{xi+i,.. . ,a;j+(„/2_i),Xi+„/2,Xj_(„/2_i), . . . (26) 

where Xi = {xi{t - si), . . .,Xi{t - s^)), Xj = Xj{t - tj) for j ^ i, f : ^ R, g : M""^ M 
are smooth functions and Tj,se € [0,r] for all j ^ i and s = l,...,m. Here, / is called the 
internal dynamics and g is the coupling function. In the papers described above, these systems 
are D^-equivariant. 

3.1.1 Characterization of delayed D„ networks 

We introduce a more general notation for delayed symmetrically coupled cell systems inspired by 
recent work on (non- necessarily symmetric) coupled cell systems of ordinary differential equations, 
see for instance [23]. Suppose that each cell in the system has phase space M*^. We generalize 
systems 1^ and 1^ to 

^iit) = fiXi,Xi+i, . . . ,Xj+(„-l)/2,^i-(n-l)/2> • • • i = 1,. . . ,n (27) 

and 

Xi{t) = f{Xi,Xi+i, . . . ,Xi^(^n/2-l)^^i+n/2^^i-{n/2~l)j ■ ■ ■ i = 1, . . . ,n (28) 

where 

Xj = {Xj{t-Ti),...,Xj{t-T^)), 

f : (M™)" — > M'^ is a smooth function and the position of X^ corresponds to the coupling from 
cell k to cell i. We say that cells j and k have identical coupling to cell i if 

f{Xi, ...,u,...,v,...) = f{Xi, ...,v,...,u,...) 

where u and v are permuted from positions j and k. We rewrite systems (I27p and (I28p as 

X = F{X) (29) 

where X = (Xi,...,X„)* 

^ = • • • >^i+(n-l)/2;-^i-(n-l)/2) • • • 

X = . . . ,Xj+(„/2-l);-''^i+n/2)-'^i-(ra/2-l)) • • • j-'^i-l) 

and F : (M™)" R"'' has z*'' component given by the formulas above for Xi(t) (for either n odd 
or even). 
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Consider the group D„, with generators p and k acting on M'^'" as follows: 

p.{Xi, . . . , Xn) = {Xn, Xi, X2, . . . , Xn~l) 

K.{Xi,. . . ,Xn) = {Xi,Xn, ■ ■ ■ ,Xn+2-j, ■ ■ ■ , -'^(n+l)/2 7 -'^{n+3)/2 ; • • • . . . , X2) if n is odd 

K.{Xi, . . . , Xn) = {Xn, . . . , Xn+i~j, ■ ■ ■ , ; ^n/2+1 ) ■ ■ ■ ,Xj, . . . , Xi) if n is even. 

(30) 

These networks can be represented by a graph as shown in Figure [3. 1.11 In this case, we have an 
eight-cell network with bidirectional nearest and next-nearest neighbour coupling with delays ti 
and T2 respectively. 




Figure 1: Representation of a Dg-symmetric coupled cell system with nearest neighbour with 
delay ri and second nearest neighbour coupling with delay T2. 

Without loss of generality we assume that the networks are transitive. That is, all cells in the 
network can be reached from any other cell via the coupling arrows. 

We now characterize the connections in the network so that the delay-differential system is 
D„-symmetric. We think of each cell in the network has having [(n — l)/2] neighbours on each 
side and an opposite cell if n is even. Graphically, it is clear that a n-cell network is D„ symmetric 
if for all cells in the network, all connections to and from the j'*'* neighbour on each side (or the 
opposite cell if n is even) are all the same; that is, the coupling term and its delay must be the 
same for all those connections. This idea is formalized in the next result. 

Proposition 3.3 A transitive network of n coupled identical cells, with delays, is Tin-e-quivariant 
if and only if it satisfies the conditions below. 

(i) Suppose that cell 1 receives an input from cell j with delay e € [0, r], then every cell i in the 
network (i = 2, . . . ,n) receives an input from cell {i + j — 1) mod n with delay e identical to 
the one received by cell 1. 

(a) For every connection in part (i), there is an identical connection from cell i to cell {i + j — 
1) mod n with delay e. 
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From Proposition 13.31 a transitive D„ symmetric network with an even number of cells must have 
nearest neighbour coupling. 

Proof: We begin by looking at p-equivariance. Denote by [w]i the i^^ row of vector w. Then, 

[pF{X)]i = f{X,.i,X„...,Xi_2). 
and since pX = Xi, . . . , Xn~i) we have 

[F{pX)], = f{Xi^i,X„...,X,_2). 

Thus, p-equivariance holds automatically by the structure of the equations. 

Consider the case n odd (the case n even has a similar proof and is omitted). Suppose that 
cell 1 receives an input from cell j. We look at the system of equation (I27p and focus on the 
possible coupling from cell {i + j — 1) mod n to cell i. Moreover, consider the possible connection 
from cell i to cell {i— j + 1) mod n. Note that the connections from (i + j — 1) mod n to i and from 
i to {i — j + 1) mod n are obtained by taking the index and subtracting j — 1. Finally, consider the 
possible connection from cell {i — j + 1) mod n to cell i. We now show that F{X) is K-equivariant 
(and so D„-equivariant) if and only if the connections defined above are identical. We see that 



/ 



kF{X) = K 



f{Xi , . . . ,Xj, . . . , Xn+2-j , 



\ 



f{Xi, . . . , . . . , . . .) 

f{Xn+2-i, ■ ■ ■ , ^n+2~i+j-l, ■ ■ ■ , Xn+2-i-j+l, ■ ■ •) 



v 



f{Xi , . . . ,Xj, . . . , Xn+2-j , 



\ 



fi^n+2-i, ■ ■ ■ ■,Xn+2-i+j-l, ■ ■ ■ i i • • •) 



/(Xj, . . . , . . . , . . .) 



and 



F{kX) 



I fi^l ) • • • ) ^n+2-j , ■ ■ ■ . . .) 

f{Xn+2-ii • • • , -'^n+2-(j+j-l)5 • • • ) -^n+2-(i-i+l) ) • • •) 
f{Xi, . . . , . . . , . . .) 



\ 
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We now show that parts (i) and (ii) imphes K-equivariance. If part (i) holds, the couphng from 
cell {i + j — 1) mod n to cell i and the coupling from cell i to cell (i — j + 1) mod n are identical. 
Then, by part (ii), the coupling from cell {i — j + l)mod n to cell i is identical to the coupling 
from cell i to cell {i — j + 1) mod n. Therefore, the coupling from cells (i + j — 1) mod n and 
(i — j + 1) mod n to i are identical. Since the dynamics of all cells is given by the same function 
/, this is true for alH = 1, . . . , n. Thus, F is K-equi variant. 

Suppose now that F is K-equi variant. Equality of both sides of the equivariance condition 
implies that for alH = 1, . . . , n, the couplings from cells {i + j — 1) mod n and {i — j + 1) mod n to 
i are identical. Since the dynamics of all cells is given by the same function /, the coupling from 
cell {i+ j — 1) mod n to cell i guarantees an identical coupling from cell i to cell {i — j + 1) mod n 
and this proves (i). But, the coupling from cell {i — j + l)mod n to i is therefore identical to 
the coupling from cell i to cell {i — j + 1) mod n. Hence there is an identical two-way coupling 
between cells i and (i — j + 1) mod n which proves (ii). I 

3.1.2 General form of the characteristic equation 

We now focus our attention on delay-coupled cell systems where each cell is one-dimensional, 
that is k = 1. We split the linear and nonlinear parts of systems l\27\ [28|) and write the result in 
abstract form: 



where Xt G C([-r, 0], M"), L : C([— r, 0], M") ^ M" is a bounded linear map and H is a nonlinear 
mapping. Thus, L is D„-equivariant, ri(9) is a n x n D„-equivariant matrix of bounded variation 
and 



(1) for all j = 1, . . . , n: r]jj{6) = p{6) for some function p, 

(2) for all i, k with i / k then r]ki{9) = ■r]{2+n~k)i{^) = %i(6')- 

Proof: We use Proposition 13.31 to obtain information on rj. By part (ii), the matrix r]{0) is 
symmetric. From the structure of (j27|28|) . we deduce that for all j = 1, . . . ,n, 'rijj[9) = p{9) for 
some function p{0). We denote by r}ji{6) the element of rj corresponding to the coupling from cell j 
to i. Consider riji{9) then there is an identical connection from cell 1 to cell 2+n—i by part (i) and 
so r]ji{9) = ^i(2+n-j)(^)- By part (ii), the connection from cell 2 + n — j to cell 1 is identical to its 
reciprocal and so: r]ji{0) = r/(2+„_j)i(6'). By part (i), we then have r]ki{0) = 'n(^2+n-k)i{G) = VkiiO) 
since the connections to cell i are identical to the connections to cell 1. ■ 



X = LXt + H{Xt) 




Proposition 3.4 The matrix rj{6) is symmetric (r]{6) = rj{6)'^) with the properties: 



Remark 3.5 This result can be obtained for higher-dimensional cells with a proof essentially 
similar to this one, but with a more cumbersome notation. We decided to restrict ourselves to 
the one-dimensional case as this is the one which we study in details in what follows. 
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For the purpose of diagonalizing the hnear equation, it is preferable to use C", the com- 
plexification of M". The decomposition of C" by the action of D„, = (p, k) yields an isotypic 
decomposition of C", see Golubitsky et. al. |13j . Let (" = e^'^*/" and define 

for j = 0, . . . ,n — 1. The subspaces Vj are irreducible representations of D„. Recall that all 
complex irreducible representations of D„ are one-dimensional. Then the isotypic decomposition 
is 

= e Fi e . . . e K-i 

where Vq is the trivial representation and for n even Vn/2 is the alternating representation both 
of which are real one-dimensional irreducible representations of D„. 
Using Proposition 13.41 and by invariance of the subspace Vj we have 



Since ('^-('^-i) = ^^-i, then for j = 0, . . . , n - 1 

n 

k=2 

n 

= p{9) + J2'^cos{2Tr{k-l)j/n)r]ki{9) 



(31) 



n odd 



k=2 

(n+l)/2 

p{e)+ 4cos(27r(A:-l)j7n)r?fciW 
k=2 

n/2 

p{9) + 2(-iyr/(i+„/2)i(0) + E4cos(2^(A: - l)j /n)rjki{e) n even 

k=2 

since ry^i = ??{2+n-fc)i aiid cos(27r(/c — \)i jn) = cos(27r((2 + n — A:) — l)j/n). Note that Aj{6) = 
An-j{9) for j = l,...,[n/2]. The block diagonalization of r] is given by the terms Aj{9) for 
J = 0, . . . , n — 1. Hence, we have 

A(A) = Mn - dri{e)e''' = XI^ - j"" diag(dAo(0)e^^ . . . ,dAn-i{e)e^'). 

Let Aj(A) = A - dAj{e)e^'^, then 

A(A) = diag(Ao(A),...,A„_i(A)). 
Therefore, the characteristic equation has, for n even, the decomposition 

(n-l)/2 

det A(A) = det Ao(A) det A„/2(A) JJ [det l^j{X)f = (32) 
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and for n odd 

(n-l)/2 

det A(A) = det Ao(A) [det Aj(A)]2 = 0. 

3.1.3 Application of Main Theorems to the D„ case 



(33) 



It is straightforward that Theorem 12.11 can be applied to any of the factors of the characteristic 
equations ([32]) and (f33l) . 

To apply Theorem 12.41 in the D„ case, we need to verify that the coefficients Irj^ in the factors 
of the characteristic equation are nonzero and that the nondegeneracy condition detX^ 7^ is 
satisfied. In fact, as it is shown in Section [H detiB 7^ if and only if the matrix 



B :-- 



1 



"1 + A*r-1 



is nonsingular where £j is the number of imaginary eigenvalues satisfying the j*^ term of the 
product of the characteristic equation (f2T]) and fij = Yli=i^i where fir = n and fj,o := 0. Note 
that row j of B contains coefficients belonging to the j^^ factor of the characteristic equation ()2ip . 
The cases n even and n odd differ significantly and we focus initially on the case n odd. 



3.2 Critical eigenvalues of D„ coupled cell systems 

In this section, we apply Theorem 12.41 to D„-symmetric coupled cell system depending on an 
arbitrary number of finite delays. 



3.2.1 n odd 

The characteristic equation is 



We can write 
where 



(n-l)/2 

det A(A) = det Ao(A) [det Aj(A)]2 = 0. 

Aj{X) = X-F{X)-Gj{X) 



F{X) = Y,aie 



Xn 



i=l 



are the terms coming from the internal dynamics of each cell and 

(n+l)/2 r ,^ .\ -, mfc 



k=2 



4cos I — — 



n 



t=i 



(34) 
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are the contributions from the couphng where ruk is the number of delayed terms in the connection 
from ceh A; to 1 and are the respective couphng coefficients. 

Example 3.6 As an example, consider a delay-coupled Ds-symmetric cell. Let Us{9) = if 
9 = [— r, —s] and Us{0) = 1 for G (~s,0] where r > s for all delays s and suppose 



p{0) V2i{0) mm mm v^m 

mm p{0) mm mm mm 

mm mm p{0) mm mm) 

mm mm mm p{0) mm 

mm mm mm mm pi.0) 



where 



P{G) = ^a.iUr,{e), ?72i (6') = ^ a^u,2 (6*) and r?3i (61) = ^ a^u^z {6) 



i=l 



with the conditions ^4i(^) = mii^) mi{^) — V2i{0) given by Proposition 13.41 part (2). Then, 



and 



k=2 



i=l 



4 COS 



n 



t=l 



where m2 = 3 and = 2. 

Thus, all coefficients bj of Is are nonzero and it is convenient to set 6{ to be the coefficient 
of oi; that is 6{ = 1 for j = 0, 1, 2, . . . , (n + l)/2, and we keep this convention for the remainder 
of the paper. 

We suppose that the characteristic equation A (A) = has purely imaginary roots coming 
from all factors, then for r = (n — l)/2 we have 



bl 



L "l Ol+Mi 



"1 + A*r-1 



-^1 + A*r-1 



and we assign the coefficients as follows. We suppose that the first row corresponds to the 

factor for the trivial representation which means that 

^'i+M. =4' J = l,2,...,(n-l)/2. 
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Then, we set the remaining coefficients of each row to be equal to 



4 cos 



27T{k-l)j- 



n 



for 



/c = 2, 3, . . . , (n + l)/2 where row j + 1 has the coefficients of Aj for j = 1, 2, . . . , (n — l)/2. This 
leads to the matrix 




1 4 cos 



[n-l)TT 



n 



4cos ( — I 



n J 

4 cos I — 

n / 



4 cos 



2(n- IW 



n 



4 cos 



4 cos 



4 

(n — 3)7r 
n 

2(n - 3)7r 
n 



4 cos 



(n — 3)(n — l)7r 
2n 



4 cos 



4 cos 



4 cos 



4 

'n — l)7r 
n 

2(n - l)7r 
n 



2n 



(35) 



Let ?! < ?2 < • • • < be a set of indices chosen from {0, . . . , (n — 1) /2} defining a combination 
of factors from the characteristic equation ()33p . We now construct the s x s matrix B by removing 
rows and columns of ([35]) not in the set {«i,«2, • • • j^s}- Suppose that ii,. . . ^ig are chosen from 
{1, . . . , (n — l)/2} then the matrix B is symmetric = B) and has the form 



B = 



4 cos 



4 cos 



27riii2 



4 cos 



4 cos 



27ri2ii 



/27ris_iii\ f 2'Kisii 

■ ■ ■ 4 cos 4 cos 



4 cos 



27ris_ii2 



4 cos 



2nisi2 



(36) 



4 cos 



4 cos 



27riiis_i 



27rziii 



4 cos 



4 cos 



27ri2is-i 



27ri2is 



4 cos 



4 cos 



/ 27^^g_l^, 



4 cos 



4 cos 



/ 27risis_i 
2ne 



In the other case, ii = and the matrix is of the form 
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4 cos 



27^^^ 



1 4 cos 



/ 2'Ki2is-i 



1 4 cos 



27ri2is 



We can now state our result. 



4 cos 



2Tris-ii2 - 1 



4 cos 



27rjs_iZs-i 



4 cos 



V n 



4 cos 



st2 



4 cos 



/ 2'iTisis 



4 cos 



V 

2^1? 



(37) 



Theorem 3.7 Consider a linear Yin-symmetric coupled cell system with n odd depending on k 
delays ti, . . . ,Tk and let ii < i2 < ■ ■ ■ < is be indices chosen from {0, . . . , (n — l)/2} defining a 
combination of factors from the characteristic equation \33\) . We assume that the matrix B given 
by ^36\) or (37\) is nonsingular. Suppose 



1 12 



are positive and linearly independent over the rationals, where li^ + ■ ■ 
exists Ti > 0, . . . , Tfc > and real coefficients Oj such that for all m = 1, 



+ 4 



k. Then there 



A.„(A) = 



has roots iujj^ for £ = 1, . . . , li^ . 
Proof: Since n is odd, the coefficients 



6{ = 4 cos 



/27r(fc-l)j 



n 



are nonzero for all A; = 2, . . . , (n+l)/2 and j = 0, . . . , (n — 1)/2. Because B is assumed nonsingular, 
Theorem 12.41 applies and the result is obtained. 1. 

The condition that B is nonsingular does not always hold as we show in the case s = 2. 
Consider the matrix (j37p with n = 9 so that i2 € {1,2,3,4}. Choosing Z2 = 3 we have the 
singular matrix 

\ 1 4 
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We now show that B is nonsingular in the case s = 2 if the matrix is given by ()36|) ; that is 



B 



4 cos 



4 cos 



27ri1 



n 



2711211 



n 



We compute 

det S = 16 



n 



27ri?\ (2TTil 
cos I cos ' 



n 



4 cos 



4 cos 



2-Kiii2\ \ 



cos 



n 



2mi 



n 



2TTiii2 



n 



27r(i? + i|) , 
cos I — ^-t — 1 + cos 



n 



2^(if 



n 



cos 



47riii2 



n 



1 



We do a few cases exphcitly. First, the case n = 3 is not relevant since ii < 12, (n — l)/2 = 1 and 
ii 7^ 0. We do the case n = 5 where we must have h = 1 and 12 = 2. This imphes that if + i2 = 5 
and so 

/47r\ /Svr' 

— cos ( 

5 



det B = 8 



cos — COS 

V 5 J 



We now turn to the general case and show that the determinant cannot vanish. Because the three 
cosines are projections of n^^ roots of unity on the real axis for n odd, then 



cos 



2Tr{il+il) 



n 



+ cos 



27r(ii 



n 



'4vriii2\ , , 
cos I ) 7^ 1. 



So if the determinant is to vanish, one of the cosines must be equal to 1. Since ii < 12 there is 
only one option and we must have if + = n. Thus, i\ = n — 12 and 



cos 



2^(if - il) 



n 



cos 



27r(n - 2il) 



n 



If 



cos 



AttH 



n 



cos 



47121^2 



n 



cos 



0, 



47ri| 
n 



this would impliy ii = 12, but we know that ii < i2 and so deti3 cannot vanish. We summarize 
this result in the next theorem. 



Theorem 3.8 Consider a linear Yin- symmetric coupled cell system with n odd depending on k 
delays ti, . . . ,Tk and let ii < i2 be indices chosen from {1, . . . , (n — l)/2} defining a combination 
of factors from the characteristic equation ^33\). Suppose 



ujI, 



1 2 
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are positive and linearly independent over the rationals, where + = ^- Then there exists 
Ti > 0, . . . ,Tk > and real coefficients ai G R, . . . , Op € M such that for m = 1 and m = 2, 

A.„(A) = 

has roots iujj^ for £ = 1, . . . , . 



3.3 n even 

For D„-symmetric systems with n even, Theorem 12.41 cannot be apphed directly because the 
condition 5], 7^ for ah indices j, k is not always satisfied. For instance, if n = 4 then 

Aj{e)=p{e) + 2{-iy7]3i{9) + Acosi27rj/4)r]2ii9) 

and we have 

Ai{e)=A3{9)=p{e)-2r]3i{e). 

Hence, = = 0. 



4 Proof of Theorem 12.41 



Before we present the proof of Theorem 12. 41 we describe in the next lemma the form of the matrix 
Zb which appears in the proof and compute its determinant. 

Lemma 4.1 Let £1, . . . ,ir be positive integers and define fij = X]i=i where fir = n and fiQ := 0. 

Consider the n x n matrix 

1b := [Ai---Aj---Arf 



where 



A, 



6^ 



is a ij X n matrix and all elements are nonzero. Then, 



-tp tp 

%+1 



detTs = ±J] 



where 



B :-- 



bi bi^^^ 



s=2 



■ b\ 



detB 



bi 



bi. 
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Proof: Substitute row k, denoted by Rk, of 



A. 



for k 



A, 



K 









L by Rk 



'^2+,! J. 



bL bi 



-bi 



■■ b^ hp -., 
Ri. The matrix Aj becomes 



-b^ b^ 



















bi 











-2y 



'2b 
















We compute the determinant of Zb by cofactor expansion starting with row 2 of Aj which contains 
a unique nonzero element —21j'^-. Denote by Cij the (z, j)-cofactor matrix. The row 2 + /ij-i of 
C(2+£j_i,/ij) has a unique nonzero element — 26_i+^^ and we perform a cofactor expansion along 
this row. The row 2+/ij„i of this new cofactor matrix also has a unique nonzero element —2ljP_2_^^. 
and we proceed with the same process removing successively columns 3 + /Xj-i to fij (and the 
appropriate rows) until the cofactor matrix has only two rows corresponding to the original Aj 
matrix and the second row has the unique nonzero element — 262+/^j_i which is used to perform 
a cofactor expansion. Performing this process successively on each matrix Aj for j = 1, . . . ,r, 
leaves as a cofactor matrix the r x r matrix B defined in the statement. The formula in the lemma 
is written using /ij = /ij-i + ij and so the lemma is proved. ■ 

We are now ready to prove our Theorem [2 



Proof of Theorem 12.41 A necessary and sufficient condition for the conclusion of the theorem 
to hold is that the following algebraic system of 2n equations has a solution in the 2n unknowns 
(Ti,r2, . . . ,r„,ai,a2, . . . ,a„): 



k=l 
n 

J2 



2g-jajlTfeN 



fc=l 



k=l 



lUJo 



ILOn 



1 £9 



(38) 
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k=l 
n 



-iwL k = 1, . . . , J 



fc=l 



(39) 



1 ■ ■ ■ i^r 



E«fc(^Ie-l-^) = 
I, fc=i 

Although (j39p is in complex form, system ([39]) is equivalent to a system of 2n rea/ equations. This 
fact is taken for granted throughout the sequel, even though we continue to use complex notation. 
Let 

u; = (u;|,...,u;£\,...,a;[,...,w;'J. 
It is useful to use the following matrix notation for (j39p 



P{-T-UJ) 



A' 



lUJ 



-lOJ 



T 



(40) 



where A = (ai,...,a„), superscript T denotes transpose, and P{t;uj) = P(ri, . . . , r„; w) is the 
n X n matrix of the form 

" Pi(t;c^) 



P(r;w) 



P2iT;u;) 

Pr{T]L0) 



whose entry at block j, row i and column k is 



Note that P(t;w) = P{-t;lo). 

Recah that T = (S^)" and V = T". Consider the following mapping associated to (I40p : 



F : V X 



defined as follows. 

F($^...,$^A;a;) 



P($\...,$'^ 



(41) 



where A and lij are as previously defined and 

Pi<^\...,<l>l 



^2(^2) 

Pri^n 
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with 



Ik 



for j = 1, . . . , r, ^ = 1, . . . , £j and A: = 1, . . . , n. The definition of P uses the fohowing coordinates 
of V. Let Tj = (§1)^:' for j = 1, . . . , r, then 



where 



= and ^i = [^{k 

Let = <J>^ and * = (^i, . . . , ^r)- We use the notation 

$ = (cl>i,...,$ri) 

where 

ci>^,= (cl>i,...,$^„_,) 
so that the mapping F in (j4ip is written as F{^, ^ , A; oj). 



^■■■'"Pikf e Tj. 



Definition: If {ei,...,e£} denotes the canonical basis of vectors in R^, we define the vectors 
vi, . . . ,V£ hy vi = ei + ■ ■ ■ + ei, and for /c = 2, . . . , 

e-2 

By construction, the set {vi, . . . ,V£} is linearly independent and so the £ x £ matrix Z^, whose j^^ 
column is the vector vj, is invertible. 

We also define the £ x £ invertible matrix Uj to be the diagonal matrix whose k^^ diagonal 
element is the k^^ component of the vector vj (in particular, Ui is the identity matrix). Note also 
that = Uj, j = 1, . . . ,r. 

We define /xq := 0, fij := Y.i=i 

We use the following base point in V = V$ x Viji. For j = 1, . . . , r define be the point given 
by 

and = -f wi for i {m+iij-i,- ■ ■ In particular, ^f,. = -^ve^ 

We now evaluate ^) by computing Pj{^, ^) for j = 1, . . . , r: 



p,($,$) = i 


'b{ ■ 
^1 ■ 


■ }r> 

■ bi,_. 


IP 


iP 


■ -bk 


tP 

%+i 


■ 




_bi ■ 








■ -K, 


tP 

%+i 


■ 
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Thus, 



Pi($,*) 



:= Hb 



where Xb is invertible by assumption. In particular, -P(— — = —Hb- We define 

^ = (ai, . . . ,a„) =T^ io 

which leads to the solution: 

F($,§,l;a;) =0. 

Because the Uj are rationally independent, it follows that the components of A are all non-zero. 
Furthermore, it is easy to compute that 



U Hb 
U -Ub 



where 



is a n X n matrix with diagonal blocks of dimensions ii x £i to £r x ir- We compute also 



K = D^F{^,^,A;u;) 



K 

k 



where 



K = dmg{Ki,...,Kr) 

is a n X (n — l)n matrix where the block Kj has dimensions ij x {n — l)£j and is of the form 
kj = (^aib{ui ■■■ ai+^ b{^^._^ Ui a2+^^._i ^2 " " " «w ^4 

The matrix J is invertible and its inverse is 



—in — 1 2T — ^ 
2 -^S 



By the implicit function theorem, there exists a neighborhood A?^ of $ in V$ and a unique smooth 
function 

G:N I — > T X M" 



G : $ 



G($) = (G*($),Ga($)) 
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such that 

G($) = ($,l) 

and 

F{^,G{^);lu) =0, £ N. (42) 
Implicit differentiation of (j42p yields that 

(43) 



where * is not important for our purposes and the first component is a n x (n — l)n matrix 
composed of n — 1 block matrices A4j = 

of dimension £j x (n — where j = 1, . . . ,r — 1. Recall that /i^ = n, so that we have 




■ t/I« t/Io . . . 5 (yli, (yl, 

-'Mr 



Consequently, 

L>Gvi;($) = diag(7Wi,...,M) (44) 
is nonsingular and it follows that the mapping 

Gm:N — ^ T 

is regular at 



Consider the following n - 1 vectors in (M^^)"-! x (M^^^n-i x • • • x (M^--)"-! ~ (M" )"-i : 

Wx = (u;^0,...,0;a;^0,...,0;...;u;^0,...,0) 
= (0,w^...,0;0,u;^...,0;...,0,w^...,0) 

VF„_2 = (0,...,a;^0;0,...,w^O;...;0,...,u;^0) 
Wn-i = (0,...,0,cui;0,...,0,u;2.^^^^0^___^0,cu'^), 
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where represents the vector in the respective space M^^ ^ and we recall that 

The set {Wi, . . . , Wn-i} is linearly independent, so for any $ G V$, the set 

5$ = { $ + ^ (7j- Wj (mod V$) I < IcTjl « 1, j = 1, . . . , n } 

is a small n — 1-dimensional surface through $ in V$. We are interested in showing that for $ 
close enough to $ in iV, the image of 5$ by in is transverse to the integral curves of the 
vector field * = To show this, we consider the function 

T -.N — >R 

defined by 

r ($) = det ( DG^ ($) • W[ DG^ ($) • W2 ■ ■ ■ DG^ ($) • ) (45) 

and recalling that (W))^ = / for all j,i we compute 



T(l') = det (^DGq,{^)-W[ DG^{^)-W^ •■ 
= det (ajk) 

where j = 1, . . . ,r, k = 1, . . . ,n. The elements of the matrix {ajk) are 



"life = < 



J^iUmuY k = l+e^,...,n-l 



(46) 



for J = 2, . . . , r — 1 



k = 1, . . . , and k = fij + 1, . . . ,n — 1 



k = n. 
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and finally 



k — 1, . . . , fij — 1 



k = n. 



where we recall that fir = n. Note that the elements of the last column are rewritten as to lead 
to the significant simplification of the determinant to the following form: 



r($) = 

where 



and for j 



detidmgiUluf, . . . Ml'\Ul)) detZ^ (47) 



T' — 

-i-R — 



Q2 
\Qr J 



Qj 



-6i 



h' 



2+//3_i 



is & ij X n matrix. Moreover, detX^ ^ since detX^ = a„ detXs and a„ ^ 0. Thus T{^) ^ 0. 



It follows that there is a neighborhood N' Q N in which T ^ 0. This is equivalent to saying 
that for all ^ E N' , the image of 5$ by in V* is transverse to the integral curves of the vector 
field * = a;. 

For each j = 1, . . . ,n — 1, the integral curves of the vector field = uj are dense in the torus 
T. Thus, for any e > 0, there is a Tj^e > and an sj^^ > such that the integral curve segment 

{ = TjLO (mod T) I Tj^e - Sj,e < Tj < Tj^s + Sj^g } 

is in the e-ball centered on ^Vj in T. For £ > small enough, the surface 

{ $ = (tiW, T2UJ, Tn-lUj) (mod V$) I Tj^e - Sj,£ < Tj < Tj^e + } 

is contained in N' and coincides with the surface <S$* for 

^* = (n.eo;, T2,£a;, . . . , Ti^ri-i),e^) (mod V$ ). 



33 



Thus, by our previous result, the n — 1-dimensional surface G^{S^*) is transverse to the integral 
curves oi ^ = uj in V^i. Since these integral curves are dense in V\ii, there are infinitely many 
intersections with Gii,(5<j)*) near the point ^ = G^{^). 

o ^ o 

Let vi/G be such an intersection point near Then there is a t„> such that 

° o 

'if=Tn (jJ (mod V^). 

o o o o o o 

Let $G 5$* be such that G^{^) =^. Then there are ri> 0, T2> 0, . . . ,r„_i> such that 

O o O O 

$= (ti T2 tj, . . . , r„_i w) (mod V$). 

o o o 

It follows from ()42|) that F(<I), vj/, G^(<I)); w) = 0, and by construction, this corresponds to a 
solution of (00]). ■ 



5 Conclusion 

We have shown in this paper that n nonresonant eigenvalues on the imaginary axis can be realized 
by a scalar delay-differential equation with n delays. Moreover, the same is true for any collection 
of n imaginary eigenvalues in a neighborhood of an n-tuple of nonresonant imaginary eigenvalues. 
We have also shown how these results can be applied to non-scalar delay-differential equations in 
the context of symmetric delay-differential equations where the characteristic equation decomposes 
according to the isotypic decomposition. We apply our result to delay-coupled D„ symmetric cell 
systems with n-odd. 

There are several ways of extending the main result of our paper. One question we did not 
address in this paper is if n nonresonant nonzero imaginary eigenvalues constitutes an upper bound 
for the realizability by a scalar equation with n delay. The case n = 1 is one such example since 
an easy calculation shows that we can have at most one imaginary eigenvalue on the imaginary 
axis. It is likely, but unknown, if this is also true for general n. 

One may want to study whether k zero eigenvalues in a single Jordan block and i nonres- 
onant nonzero imaginary eigenvalues can be realized in a scalar delay-differential equation with 
k + i delays. This problem may be feasible by modifying the proof of Theorem 12.11 since the 
nonresonance of the i eigenvalues is again present. However, we expect the argument used in this 
paper to breakdown for n nonzero imaginary eigenvalues with resonance. We can also study the 
same problem as in this paper but for higher dimensional delay equations. One problem would 
be to find out if n nonresonant nonzero imaginary eigenvalues can be realized by a m dimensional 
system with k delays. For instance, it is known that a pair of nonzero imaginary eigenvalues can 
be realized by a two-dimensional equation with one delay [8j. In this case n = 2,m = 2,k = l and 
so 71 = mk; is it possible to realize three nonresonant nonzero eigenvalues or does the relationship 
n = mk provide a bound to realizability in general? 

Another problem which can be studied is whether a restriction in the class of delay equation 
can change the realizability. For instance, the characteristic equation for a general two-dimensional 
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system with one-delay r is 

X^ + aX + bXe-^^ + c + de'^^ = 

while for a second-order equation with one delay in the feedback term we must set 6 = 0. We 
know in this case that two nonresonant nonzero imaginary eigenvalues can be realized by a second- 
order equation with a unique delay in the feedback term [3]. In fact, two imaginary eigenvalues 
with 1 : 2 resonance has been found in such an equation [5]. The obvious question is to see if n 
nonresonant (and resonant) nonzero imaginary eigenvalues can be realized within the class of n^^ 
order scalar equations with one delay in the feedback term. 

An extension of our main result in a direction relevant for studying bifurcations is to find 
out whether the n nonresonant nonzero imaginary eigenvalues can be realized by a scalar delay 
equation such that the remaining eigenvalues have negative real parts; that is, the multiple Hopf 
point lies at the boundary of the stability region for the equilibrium solution. 

Finally, let us mention the case of linear T-periodic equations with + 1 delays: 

N 

X = '^aj{t)x{t - Tj) (48) 

j=0 

where each aj{t) is a T-periodic n x n matrix. Hale states the following open problem: 

Is it possible to give a precise upper bound in terms of N on the number of Floquet multipliers 
of ^8\ ) that can have moduli 1 ? 

If we restrict equation (08]) to scalar equations we can pose a possibly simpler problem which is 
related to the main result of our paper: 

Is it possible to realize + 1 complex numbers e^*"^^ , • • • , e^*^^+^ with loi, . . . , io^+i positive and 
rationally independent as Floquet multipliers of the scalar equation ^8\ ) ? 

This problem is automatically solved by Theorem 12.11 if a Floquet theorem can be applied 
to (08]); that is, the Floquet exponents of the Floquet multipliers of equation (08]) are eigenvalues 
of a scalar equation ([T]) with A^ + 1 delays. Such a theorem has not been proved in general, 
however it may hold true given some conditions are imposed on (08]). 
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